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Q Based on observations collected at the SEST 15m tele- 
scope (Chile). 
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ABSTRACT 

We present 1.25 mm continuum data for a southern galaxy sample selected from the 
IRAS PSC and complete to S^o — 2 Jy. Two thirds of the galaxies have been detected 
and significant limits on the remaining objects have been set. We find, on a statistical 
basis, indications that the dust emission in these galaxies is somewhat more centrally 
concentrated than that of the optical light, possibly tracing a higher metal content in 
the inner galactic regions. This result also allows to estimate the aperture corrections 
to the millimetric data. The latter, together with IRAS photometric data, have been 
used to compare the broad-band FIR/mm spectra with a simple dust model. According 
to their far-IR/mm spectrum, the sample galaxies show a dichothomy: almost half of 
the objects, those displaying bright 25-60 fim fluxes ascribed to warm dust residing in 
starburst regions, are characterized by higher values of the bolometric (optical + FIR) 
luminosity, of the dust-to-gas mass ratio, of the dust optical depths and of the overall 
extinction. A complementary class of objects dominated by cold dust (cirrus) shows 
opposite trends. Because of the favourable observational setup, selection wavelength 
and completeness, we believe these data provide an exhaustive and unbiased view of 
dust properties in spiral galaxies. 

Key words: galaxies: ISM - galaxies: photometry - galaxies: spirals, starburst - ISM: 
dust, extinction - infrared: ISM: continuum 



1 INTRODUCTION 

Observations of diffuse dust in galaxies provide important 
information on some basic questions about their present 
structure and past history. Even minor amounts of absorb- 
ing dust play a crucial role in shaping the galaxy's broad- 
band e.m. emission. Any inferences about structural prop- 
erties, such as the global baryonic content, its spatial dis- 
tribution, and the presence of a dynamically dominant non- 
baryonic component, rely on our ability to account for the 
missing fraction of light in optical searches (Disney, Davies 
and Phillips, 1989; Choloniewski 1991; Burstein, Haynes & 
Faber 1991; Huizinga & van Albada 1992; Giovannelli et al. 
1994; Valentijn 1994; Peletier et al. 1995). Are we missing 
significant amounts of luminous matter because of the effects 
of dust extinction? How much severe are the corresponding 
selection effects? 

Dust is also an important tracer of the activity of stellar 
populations in galaxies. By cumulating the products of all 
stellar generations, diffuse dust provides an integrated view 
of the past history of the object. 



There is a tendency, as of today, to consider the inner- 
most parts of late-type galaxy discs close to optical thick- 
ness, but for the outer regions the question cannot yet be 
settled with present data. Fundamental questions of extra- 
galactic astronomy, e.g. the calibration of distance indica- 
tors, are affected by this lack of knowledge. 
Traditional approaches to investigate dust in galaxies rely on 
either (a) indirect estimates based on dust extinction effects 
in the optical and near-IR (e.g. White & Keel 1992; James 
& Puxley 1993; Block et al. 1994), or (b) direct measure- 
ments of dust emission from FIR/mm observations. Neither 
have provided conclusive results at the moment. The former 
method is based on rather model-dependent and controver- 
sial assumptions about the broad-band optical-IR spectrum 
of various galactic components. The latter method, on the 
other hand, still lacks a large enough statistical basis and 
suffers from the observational uncertainties of the mm data. 

IRAS survey data, although having provided the start- 
ing point for studies of galaxy emission at long wavelengths, 
are not enough by themselves to fully characterize the dust 
content and emission properties, because of its limited sen- 
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sitivity, spatial resolution, and spectral coverage. In partic- 
ular, dust colder than 20 K at large radial distances, whose 
existence is currently matter of so much controversial dis- 
cussions (Rowan-Robinson, 1992; Valentijn, 1994; Block et 
al., 1994), can only be sampled at A > 100/im (Chini et 
al., 1986; Thronson et al., 1987; Stark et al., 1989; Andreani 
and Franceschini, 1992), where direct observations are still 
limited to a handful of objects. 

This paper reports on a millimetric continuum survey of 
a complete sample of IRAS galaxies performed at 1.25 mm 
with the 15m Submillimeter ESO-Swedish telescope (SEST). 
A preliminary report on these observations has been given 
by Franceschini and Andreani (1995). 

The SEST telescope was chosen as it provided the best com- 
promise between detector sensitivity and spatial resolution. 
The 24" FWHM of the SEST beam at 1.25 mm and the aver- 
age optical extent of the sample galaxies match favourably, 
so that the beam-aperture corrections needed to compare 
with the IRAS survey data are not as severe as for other 
observational setups. 

Observations at such long wavelengths, made possible by the 
high sensitivity of the used bolometer system, ensure that 
all possible dust components with temperatures down to the 
fundamental limit of T ~ 3 K set by the cosmic background 
radiation are properly sampled. The reference to a complete 
flux-limited sample ensures a minimal exposure to the ef- 
fects of bias, hence allowing any kinds of statistical tests to 
be applied (determinations of IR-mm bivariate luminosity 
distributions, luminosity functions and volume emissivities: 
Franceschini and Andreani, 1996, work in progress). Finally, 
the adoption of a far-IR selected sample, rather than of an 
optical one, makes us confident that the whole phenomenol- 
ogy of dust effects in galaxies is properly sampled. 

Section 2 is devoted to the present observations and 
discusses the issue of the extent of the millimetric emission 
and the beam-aperture corrections. Section 3 reports results 
from the best-fit to the observed spectra with a dust model 
and the inferred dust properties for our sample, while the 
discussion is deferred to §4. 



(LED A), the NASA Extragalactic Database (NED) and the 
IRAS PSC (Moshir et al., 1989). Optical photometric data 
have been taken from the ESO-LV catalogue (Lauberts & 
Valentijn, 1989). 



Table 1 



2.2 The bolometer observations 

The sample has been observed at 230 GHz during various 
campaigns from 1990 to 1993 with the SEST 15m telescope 
feeding a 3 He-cooled bolometer of the MPIfR (Kreysa 1990). 
The diffraction- limited beam size is 24" (FWHM). Beam- 
switching is achieved by chopping ON-OFF the source and 
nodding the telescope, which results in a three-beam modu- 
lation on and off the source. The beam separation was set at 
70". The pointing accuracy was most of the time 2" and was 
checked each half a hour, by pointing a nearby radio-loud 
quasar. Each source was observed n x 200 seconds, with 
n depending on the expected 1.25 mm intensity (approx- 
imately evaluated by extrapolating the IRAS 100pm flux 
using a modified thermal spectrum with spectral indices of 
dust opacity between 1 and 2). 

The atmospheric transmission was monitored by fre- 
quent skydips. Uranus, Mars and Neptune were used as pri- 
mary calibrators and the pointing quasars as secondary ones. 
The overall accuracy of the measured fluxes (which we es- 
timate to be between 10 and 30%) strongly depends on the 
atmospheric conditions. A more detailed description of the 
observing procedure can be found in Andreani & Frances- 
chini (1992) and Andreani (1994). 

Table 2 reports our observed 1.25 mm fluxes (cols. 4 and 
5), together with the IRAS at 25, 60 and 100 (im fluxes (cols. 
1, 2, 3), all corrected for color and K-correction following 
Smith et al. (1987). In the lack of a detection a 3 a upper 
limit to the flux is reported. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 The sample 

The galaxy sample under investigation was selected from 
the IRAS Point Source Catalogue and is complete to the 
flux limit of S60/jto = 2 Jy. It consists of 30 galaxies with 
morphological types from SO/a through Scd, within the sky 
region 21 h < a(1950) < 5 h and -22.5° < S < -26.5°. For 
most of these galaxies PSC data also at 100, 25 and 12 \im 
were available. Distances (mostly from distance indicators, 
otherwise derived using redshift data for Ho = 75) cover 
the range from d ~ 20 to ~ 200 Mpc. Half-optical-light di- 
ameters A e are typically within A e ~ 10" to ~ 40", four 
exceptionally close galaxies having A e ~ 100". The volume 
test (< V/Vmax >= 0.45 ± 0.06) ensures the sample com- 
pleteness. 

Table 1 reports some relevant information (name, coor- 
dinates, distances, optical dimensions, blue magnitudes and 
morphological types) for the sample galaxies. The data have 
been taken from the Lyon-Meudon Extragalactic Database 



Table 2 



2.3 Beam-aperture corrections: the size of 
galaxies at A = 1.25 mm 

Corrections to the observed millimetric fluxes, to account 
for the beam aperture of the telescope being comparable 
to or even smaller than the angular size of the sources, are 
crucial to estimate the global mm emission and to compare 
it with the total fluxes provided by the IRAS survey. Such 
a difficult evaluation is one of the main reasons preventing 
so far to reach reliable conclusions about dust emission in 
galaxies. 

Up to now our knowledge on the dust distribution in 
galaxies mostly relied on IRAS data, whose basic results are 
briefly summarized in the following. 

Light profiles similar, on average, in the far-IR and opti- 
cal band were inferred from a direct comparison of blue- 
band, A = 50 and 100 [im data using the IRAS CPC 
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(Chopped Photometric Channel) instrument by Wainscoat 
et al. (1987) on three nearby edge-on spirals and by van 
Driel et al. (1995) on 55 nearby spirals. A study of a sam- 
ple of extended galaxies, partially resolved by IRAS (Rice et 
al. 1988), shows that the mean ratio of the far-IR (60 /im) 
Die, to blue-light Db isophotal diameters turns out to be 
0.98±0.25. This implies that on average galaxies have far-IR 
extensions comparable to the blue ones. A similar result is 
found by Xu & Helou (1995) for M31. They combined FIR 
and UV data of M31 and a dust heating/cooling model, to 
study the extinction in the disc of this galaxy. The over- 
all distribution of V-band optical depth in this galaxy is 
rather flat and it is very tight correlated with the HI col- 
umn density. There are IRAS objects, however, where the 
IR emission seems to be more centrally concentrated than 
that of the blue light. 

First attempts to map mm dust emission have been 
recently done by Chini & Kriigel (1993) and Chini et al. 
(1995), who have shown that the distribution of the milli- 
metric light is similar to the optical one and argue that for 
the observed objects cold dust emission is extended as far 
as the blue light. 

Two main advantages distinguish our sample: the ma- 
jority of the objects are relatively distant and have small 
angular sizes in the optical. Furthermore, the SEST antenna 
has a a relative large beam-width with respect to other milli- 
metric telescopes (JCMT and IRAM). In spite of this, aper- 
ture corrections may still turn out to be significant and de- 
serve particular attention. The solution we devised is the 
following. For some of the optically most extended sources 
(0300-23, 0335-24, 2217-24, 2233-26) we have performed 
a rough mapping by integrating on four points, one beam- 
width distant, around the galaxy center. The fluxes listed 
in Table 2 have been obtained by adding the contributions 
of each observed position. This allows at least to reduce the 
large aperture corrections for these sources. 
The galaxy 0045-25 (NGC 253) was not observed by us and 
the flux listed in Table 2 is taken from Chini et al. (1984). 
These authors performed aperture photometry with a beam 
of width of 3.9' and measured a flux of 4.4±1.7 Jy. The inner 
region (2'x3') of this object was also mapped at 1.25 mm by 
Kriigel et al. (1990). These authors estimate a total flux in 
this inner region of 1.7±0.2 Jy. However, this well-known 
object is very nearby and therefore quite extended in the 
optical (27.4'x6.8'), it is very likely that both measurements 
have underestimated the real emission. 

More generally, our adopted procedure was to test three 
different hypotheses about the geometrical distribution of 
the millimetric flux. (1) Sources are point-like objects, hence 
no beam-aperture correction is to be applied; (2) the ra- 
dial distribution of the millimetric light is exponential with 
a scale-length a mm equal to one third of the optical one: 
Qmm = ao /3 (case of dust emission more concentrated than 
the star-light) and (3) the dust emission follows that of the 
optical light (a mm = ao). 

The model surface brightness distributions corresponding to 
the three cases have been convolved with the telescope beam 
shape following Andreani and Franceschini (1992). Results 
of these aperture corrections to the observed mm fluxes are 
reported in Table 2. In columns 6 and 7 the fluxes corrected 
according to hypothesis (2) are listed, and in columns 8 and 
9 those corrected according to hypothesis (3). Note that for 



Figure 1. (a) Millimetric to far-IR flux ratios versus galaxy dis- 
tance. Top panels compare with the 100 fim flux, bottom panels 
with the 60 fim. In panels to the left the mm data have been cor- 
rected for beam-aperture assuming point-sources, central panels 
assume a mm scale-length of 1/3 of the optical one, right-hand 
side panels assume the same distribution for the mm and optical 
light. Panel (b) same as (a), but for the far-IR flux ratios versus 
galaxy diameter. 



case (2) the aperture-correction turns out to be 20% on av- 
erage of the observed mm flux (the largest correction is a 
factor 3 for 0311-25). An average correction of roughly a 
factor 2 is instead required according to hypothesis (3). 
We have tested further the three different hypotheses on the 
mm light distribution against our data as follows. Figures la 
and lb are plots of the ratios of the 1.25 mm to the 60 and 
100 fim fluxes versus galaxy distance and diameter, with 
Si. 25 corrected according to the three above mentioned hy- 
potheses. The two lines in each panel of Fig. 1 correspond to 
the ±1 a regression fits obtained with a "survival analysis" 
technique. 

We see in Fig. la that hypothesis (3) in particular im- 
plies a dependence of the average flux ratio with distance 
which is significant at the 3.5cr level for the 100 fim (and 
5.5a for the 60 fim match). This is of course not expected, 
since all flux measurements have already been K-corrected. 
We are led to conclude that the hypothesis of the presence 
of emitting dust at large radial distances appear to be incon- 
sistent with our data, as detailed in Sect. 3 below. Instead, 
the case for a higher concentration of dust emission, with 
Qmm between cv /3 and a /2, is fully supported. An even 
stronger indication in favour of the latter assumption is ob- 
served in Fig. lb, where the mm-IR flux ratios versus the 
galaxy effective diameter is shown. 
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The millimeter fluxes used hereafter were therefore corrected 
for the effect of beam-aperture according to the assumption 
a mm = ao/3. Our mm data for all the sample galaxies, as 
well as the complementary IRAS fluxes, are plotted in Figure 
2. Best-fits of a dust emission model discussed in the next 
Section also appear in the figure. 



3 THE DUST CONTENT OF GALAXIES 
3.1 The dust model 

Data at millimeter wavelengths, improving on the spectral 
dynamic range with respect to the IRAS survey data, add 
significantly to our knowledge of dust properties in infrared 
galaxies. No dust component, at whichever temperature, is 
expected to escape a combined mm-IRAS survey. 

FIR spectra of spiral galaxies are schematically inter- 
preted as due to thermal re-radiation from dust present in 
three main galactic environments: (a) dust grains heated 
by the general interstellar radiation field ("cirrus" emission 
including cold dust and small transiently-heated hot dust 
grains), (b) dust heated by hot stars in extended HII re- 
gions, and (c) dust in compact HII regions and molecular 
clouds. The "cirrus" component has been modeled follow- 
ing the recipes by Rowan-Robinson (1986, 1992; hereafter 
RR86, RR92). The model predicts the existence of nine 
grain types: amorphous carbon and silicate grains of size 
a — 0.1 /ttm, graphite and silicate of 0.01 and 0.03/im size, 
very small (a = 0.002 and 0.0005 /im) transiently-heated 
graphite grains, and finally a component of very large amor- 
phous grains with radius a = 30 /Km to account for the 
observed submillimeter properties. The mixture and grain 
properties have been optimized by RR92 to reproduce all 
basic observables (extinction law, absorption properties of 
the ISM towards various galactic regions) over the entire 1 
mm to 0.1 /itm spectral range. 

For normal grains the emission temperature is assumed 
to attain equilibrium (the thermal capacity is large enough 
with respect to the energy of the field photons), and is esti- 
mated from an energy balance equation between the emitted 
and absorbed photons (including the scattered fraction). On 
the contrary, very small grains, because of their small sizes, 
are not in thermal equilibrium with the impinging radiation 
field and their emission fluctuates according to their changes 
in temperature. Their thermal behaviour is described by a 
probability distribution function P(T), which has the ap- 
proximated analytic expression given by Draine & Anderson 
(1985): 



p( T) = { 6.680 ■ T- 2 - 75 if 2.7 < T < 80 a = 20A; 
1 ' \ 0.168 ■ T~ 2 - 75 if2.7<T< 500 a = 5A 

Their scattering and absorbing properties for A > 0.1/xm 
are those of 0.01/im grains scaled by a factor of a/0.01 (see 
below). In fact for small grains (a < 0.1 A) the extinction 
efficiency, Qext/a is independent of a. 

For the warm dust in star-forming regions we have fol- 
lowed a simplified but successful approach as illustrated by 
Xu and De Zotti (1989) and Conte (1993). The star- forming 
region is modeled as a volume uniformly filled with hot stars, 
dust and radiation, and is assumed to be optically thin at 
least for A > 20 fj,m. The same grain mixture as for the 



Figure 2. FIR/mm spectra of the sample objects. In each panel 
the name, the total FIR/mm luminosity and dust mass (estimated 
from the dust model) are shown. Solid lines are the global fitted 
spectrum, including two dust components: (1) the cirrus (dashed 
line), and the starburst component (dot-dashed line). 
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"cirrus" component is adopted, but very small grains are 
assumed to be destroyed by the intense radiation field. 

3.2 The dust emissivity 

The total emissivity due to the hot, warm and cold compo- 
nents is given respectively by: 

Yl Xc J ^ N 9r,H(ai)Q ext {ai, v)nalP{T)B u (J)dT+ 

i— 1 

4 

^^N ar , H {a i )Q £xt {a l ,v)-Ka 2 l B v {T,)+ (1) 

i=l 
7 

^ 4irN gTl H(a,i)Qext(ai, v^a^B,^) 

i=l 

N gr ,H are the number of grains per H atom, Q ex t(cn, v) are 
the extinction efficiencies, B v (Ti) is the black body spectrum 
at the equilibrium temperature and \c is the ratio between 
the intensity of the interstellar radiation field and that of 
the field in the solar neighbourhood. This latter has been 
taken by Mathis et al. (1983). The sum of the warm compo- 
nent runs over only four of the seven components, because 
for warm environments, such as those of HII regions, the 
differences in optical properties of grains with equal dimen- 
sions and different chemical compositions can be neglected. 
In this case, the averaged extinction efficiencies and total 
abundances have been put in eq.(l). The UV absorbing effi- 
ciency for these grains was taken equal to 1 and the radiation 
field was assumed uniform inside the HII regions. 

The equilibrium temperatures for all grains of both the 
"cirrus" and starburst components are uniquely determined 
by solving the balance equations with the the average radia- 
tion field intensity, Xc, in the cold "cirrus" and the intensity 
illuminating the warm dust in star-forming regions, \w, as 
free parameters. The third parameter used in the model is 
the light fraction f w at 100 fim contributed by warm dust to 
the galaxy spectrum. These three parameters fully describe 
the dust model. The fitting procedure assumes that the 1.25 
mm flux is dominated by cold dust belonging to the "cir- 
rus" component, which mostly defines \ c . The 60 and 25 
fxm data mostly define the other two parameters. 

3.3 Dust distribution in galaxy discs 

We test in this section which spatial distributions of dust in 
galaxy discs correspond to the surface brightness distribu- 
tion of the millimetric emission inferred in §2 from our ob- 
servations (though we recall that this distribution has been 
indirectly estimated using a statistical argument, and then 
may be prone to unaccounted systematic effects). We ex- 
pect, in general, that because of the scaling with radius of 
the interstellar radiation field and hence of the dust emis- 
sivity, the mm light should be more centrally concentrated 
than the optical. 



Figure 3. Integrated galactic light as a function of the galactic 
radial distance. The solid curve corresponds to the starlight, the 
dashed curve to the dust emission found by assuming the dust 
distribution follows that of the stars. The dotted curve represents 
the integrated dust emission with a dust distribution more con- 
centrated that that of stars (the scale-length in this case is one 
third of the optical). 

We have modelled the galactic interstellar radiation field 
as an exponential scaling with the radial distance in the 
disc I = Io ■ exp(—r/rd) (with rd~3.5 kpc, e.g. Reid 1994), 
i.e. closely following the distribution of stars. We have then 
computed the dust emission according to our adopted dust 
model, assuming (a) that the dust distribution is more cen- 
trally concentrated (with scale-length of one third of that 
of stars) and (b) alternatively that it follows the starlight 
distribution. 

The integrated surface brightness distributions resulting 
from this calculation are shown in Figure 3, the solid curve 
referring to the interstellar radiation field, the dashed one 
to a mm emission resulting from a dust distribution equal 
to that of stars, the dotted curve corresponding to the case 
of a dust scale-length of one third of that of stars. The cor- 
responding ratios of the optical scale- length, a , to the mm- 
light scale-length, a mm , are 1.3 in the first case and 3.1 for 
the second. Obviously, for a centrally concentrated dust dis- 
tribution, the key factor is a mm rather than a . This implies 
that the spatial distribution of the mm light that we have 
inferred roughly coincides with the real underlying dust dis- 
tribution. 

3.4 Best-fit dust models to the observed spectra 

The best-fit spectra for the sample galaxies are shown in 
Fig. 2. The adopted dust model, for a suitable choice of the 
parameters, reasonably reproduces the observed far-IR/mm 
broad-band spectra for all objects, a result which confirms 
and extends previous findings by RR92. The dashed line 
in Fig. 2 corresponds to the "cirrus" component, the dot- 
dashed ones to the starburst component. 

For 16 galaxies in our sample the contribution of warm 
dust at 100 jxm appears to be negligible, and sometimes not 
even required to improved the fit with a pure "cirrus" com- 
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ponent. We call these " cirrus" -dominated galaxies. On the 
other hand, 14 objects require high values of the warm dust 
fraction (f w > 0.3 0.4), and their spectrum at A < 100 fim 
is dominated by the starburst component (the starburst- 
dominated galaxies). We then suggest that far-IR to mm 
data and a reliable dust model provide a straightforward 
way to classify galaxies according to their star formation ac- 
tivity. This allows a complementary classification tool with 
respect to the usual approach based on optical-UV spectro- 
photometry. 

In each panel of Fig. 2 best-fit values for the total dust 
mass and the FIR luminosities are also reported. The lat- 
ter have been computed by integrating the total FIR/mm 
spectrum. Dust masses for both cold and warm dust have 
been estimated following RR86's recipes for the mean den- 
sity of grain material, while, for a given spectral fit, the total 
number of grains is set by the far-IR flux normalization. The 
mass of diffuse dust in the interstellar medium is found to be 
mostly contributed by the cold "cirrus" component (97.5% 
on average), with only a minor fraction of the mass resident 
in the starburst component. The dust fraction involved in 
the latter ranges from typically several percent for starburst 
galaxies to one percent in the inactive population. Figure 
4 shows the average FIR/mm spectrum for 30 galaxies in 
our sample fitted with a two-dust component and the cor- 
responding best-fit parameters of the dust model are listed 
in Table 3 (the reduced xt is 0.68). Note that the starburst 
component contributes only a couple of percent on average 
of the total dust mass. 



Table 3 



The quoted dust masses should be taken as lower lim- 
its: some dust may be hidden at low enough temperatures, 
given the strong dependence of dust emissivity on tempera- 
ture. We argue here that only very sensitive mm observations 
(like those presented here) can tackle it successfully. In fact, 
though the results from these fits are by no means unique, 
there is no room for large amounts of dust mass in the frame- 
work of the adopted model. The data are compatible with 
models having an additional very cold dust component il- 
luminated by a very dim radiation field (x c = 0.5x©) and 
with a temperature distribution ranging from T~4 to 18 K. 
In this case this very cold component would bear 76.3 % of 
the total mass and the total average dust mass would be 
increased by 30%. 

The largest amount of dust still allowed by our observed 
dust spectra may be computed assuming that all dust com- 
ponents are set to the fundamental temperature of 3 K (the 
corresponding radiation field would be \c = O.Olx©)- The 
resulting mass of this very cold dust would include 90% of 
the total, which would be twice larger than our best-fit value. 
Then we estimate that our dust masses are correct within a 
factor of two, for the adopted dust model. 



Figure 4. The average far-IR millimeter broad-band spectrum 
for the median galaxy in our sample. A two component model, 
including the cirrus (dashed line) and starburst (dot-dashed line) 
contributions with best-fit parameters \c = 11-X0: X™ = 120xo 
(normalized to solar neighborhood) The flux normalization at 100 
fim of the starburst component is f w = 0.3. The physical param- 
eters of dust grains are detailed in Table 3. The median value of 
the 100 fim flux for our sample galaxies is Sioo = 7.5 Jy. 



3.5 Effects of dust extinction on the overall 
spectrum 

An overall view of dust effects in our IR galaxy sample is 
summarized in Figure 5 in terms of extinction versus dust 
optical-depth. The observed B-band optical depth tb has 
been estimated from the total dust mass divided by the pro- 
jected area within which we estimate to be resident the cold 
dust, i.e. the inner third of the optical radius (R25), prop- 
erly scaled with the axial ratio squared (b/a) 2 to account 
for the disk inclination, tb, which measures the amount of 
dust available to absorb the optical light, is compared in Fig. 
5 with the quantity Ab measuring the overall actual effect 
of extinction. The average B extinction, Ab, has been esti- 
mated from the logarithmic ratio of the bolometric optical- 
UV luminosity from 0.1 to 10 fim (Lo) to the bolometric 
far-IR light (Lfir, 10 to 1000 fim) (two objects are missing 
in Figure 5 because optical photometry is not available). 

We see that the two galaxy classes, that we have de- 
fined are significantly segregated over this plane, the in- 
active "cirrus" -dominated objects being confined to lower 
values of dust optical depth (tb < 2) and low extinction 
{Ab < 1). The active star- forming galaxies, on the contrary, 
are spread over much larger values in both axes. The pre- 
dicted dependences for a screen, a slab, and a sandwich (with 
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Figure 5. Observed optical depth of dust Tg versus the overall 
extinction Ag. Tg is averaged within one third of the R25 radius. 
Ag is estimated by comparing the bolometric outputs in the op- 
tical and far-IR. Open and filled squares refer the the inactive 
cirrus-dominated and to the starbursting objects. The position 
marked by the symbol (0) corresponds to our Galaxy. 



Figure 6. Best-fit of the dust model to spectral intensity data on 
our Galaxy (Pugct 1989; Page ct al. 1990; Wright et al. 1991). The 
data are fitted by a pure cirrus emission (no starburst emission), 
with Xc = 8X0- The total dust mass is M d = 10 7 Mq. The 
dashed line corresponds to the best-fit by Wright et al. (1991). 



zero scale-height of dust) model are shown for comparison 
(see Disney et al., 1989). Appreciable amounts of dust and 
extinction seem to characterize only a population of star- 
bursting galaxies, with typical values of the effective optical 
depth Tg of 1 to 20 in the inner one third of the optical ra- 
dius. However, only for a minority of these starbursts (4/13) 
extinction values significantly higher than 1 are indicated. 

The solar symbol (0) in the plot of Fig. 5 marks the 
position of the Galaxy. We have estimated it by applying 
the model of §3.1 and §3.2 to the Galactic integrated spec- 
tral emissivity, as shown in Figure 6. The Galaxy spectrum 
is well fitted by a pure cirrus emission with Xc/x© — 7. The 
estimated total dust mass turns out to be 10 7 Mq while 
the optical depth r ~ 0.3. These values are somewhat lower 
than estimated by Sodroski et al. (1994), because of our dif- 
ferent fit to the long-wavelength part of the spectrum (see 
comparison Fig. 6). The position of the Galaxy within the 
region occupied by the inactive galaxy population is consis- 
tent with our finding that the Galaxy's far-IR spectrum is 
dominated by cold cirrus emission. 

3.6 The dust-to-gas mass ratio 

We report in Figure 7 a plot of the gas (HI + molecular, 
from Andreani et al., 1995) to dust ratio versus bolometric 
luminosity (L bo i — Lo +Lfir). We first notice that the me- 
dian Lboi ~ 4.5 10 10 Lq for the inactive objects is a factor 
2.6 lower than the corresponding value for starbursts, consis- 
tent with our inference that the latter are characterized by 
an enhanced star-formation activity, hence higher luminosi- 
ties. Second, the median gas/dust ratios for the two classes 
also differ on average (< M g /Md >^ 400 for the starbursts 
and ~ 1000 for the inactive population). Finally, there is an 
apparent trend of < M g jM d > to decrease with L bo i from 
values of 3000 to less than 400 for the inactive population. 
The gas-to-dust ratio is particularly large (>3000) for 



Figure 7. The gas to dust mass ratio versus bolometric opti- 
cal/mm luminosity (symbols as in Fig. 5). Mass of gas (including 
both HI and molecular gas) estimated from CO observations (see 
Andreani ct al., 1995). 

three quite extended and close-by objects: 0300-23, 0045-25 
(NGC 253) and 0128-22. The most likely explanation is that 
we have missed part of the dust emission, or underestimated 
the aperture correction to be applied (see §2.3), for 0300-23 
and 0128-22, and it is probably even more so for the obser- 
vations of NGC 253 by Chini et al. (1984) and by Kriigel 
et al. (1990), whose photometry refers only to the inner re- 
gion. In principle, it could be argued that these three objects 
are truly metal poor, but this is unlikely since CO emission 
were detected in all of them (Andreani, Casoli & Gerin 1995; 
Devereux & Young, 1990). 

The position of the Galaxy in Figure 7 has been found 
by scaling the value of the ratio < M g /Md > of Sodroski 
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ct al. (1994) to be consistent with our lower estimated Ma- 
lt is confirmed that the Galaxy tends to have a somewhat 
low gas-dust ratio, but the difference with respect to the 
median value of inactive galaxies of the same luminosity is 
now a factor 2, only marginally significant. 



4 DISCUSSION 

A fruitful combination of good sensitivity at long wave- 
lengths, a large beam aperture and small angular size of 
the target objects, have allowed us to explore the dust con- 
tent in galaxies down to the coldest possible temperatures 
and over a significant fraction of the optical galaxy extent. 

Although our inferences on the millimetric size of galax- 
ies have only a statistical sense, our results in Fig. 1 ap- 
pear inconsistent with the assumptions that cold dust emis- 
sion has a scale-length comparable to or larger than that of 
starlight. They rather suggest that the mm scale-length is 
somewhat smaller than the optical one, though is not possi- 
ble to precisely quantify the difference until large bolometer 
arrays will be available. 

From a detailed modelling of dust emission from a galaxy 
in which the dust is illuminated by an exponentially dis- 
tributed radiation field, we conclude that our inferred dif- 
ference between optical and mm lights corresponds to a dust 
distribution which is more centrally concentrated than the 
starlight, with a typical scale-length of one third of that of 
stars. 

This effect may bear some relationship with the observed 
increase in metallicity towards the inner regions of galaxies 
(e.g. Issa et al. 1990; Sodroski et al., 1994; Carollo, Danziger, 
Buson, 1995) and with the idea of an enhanced past stellar 
activity there. A similar effect is seen in the Galaxy and 
in M31, where a dust-to-gas ratio increasing towards the 
inner Galactic regions is observed (Sodroski et al. 1995; Xu 
& Helou 1995). This fact is related by Sodroski et al. to 
the metallicity gradient as a function of the galactocentric 
distance (Panagia & Tosi, 1981; Shaver et al. 1983; Wilson 
& Rood 1994) observed in the Galaxy. 

A published model (RR86, RR92) reproduces quite well 
the observed broad-band spectra (see Figures 2 and 4). Ac- 
cording to their FIR/mm spectrum, a simple color criterion, 
supported by model predictions, allows us to classify these 
IR galaxies in two different groups: an inactive galaxy popu- 
lation, the we called cirrus dominated objects, and an active 
star-forming population, called starbursts. 
The average galaxy FIR/mm spectrum can be fitted by a 
dust model including contributions from a cold and warm 
dust components, with the cold one bearing most of the dust 
mass. We find for the median galaxy in our sample a value 
of 1.9 10 7 Mq in dust, only 2 percent of which is warm dust 
associated to star-forming regions. 

We have tested the presence of very cold dust, still in the 
frame of the adopted dust model, by putting all dust compo- 
nents to the fundamental temperature of 3K (illuminating 
it with a very low radiation- field intensity). This, probably 
extreme, condition would increase the average dust mass to 
3.5 10 7 Mq, 90% of which in very cold dust. We argue here 
that in the frame of the adopted dust model there is no 
large room for hidden cold dust and the estimation of the 
dust mass is quite robust. 



It must be noted that estimation of the dust masses 
based on IRAS data only would produce values smaller by 
one order of magnitude. This was already noted by other 
authors. For instance, even indirect methods used to infer 
the amount of dust in galaxy disks, such as those based on 
the observed B-K and V-K colours, together with population 
synthesis models (Block et al. 1994), or the combination of 
FIR and UV data with a dust heating/cooling model (Xu 
and Helou, 1995), agree in finding dust amount larger than 
that inferred on the basis of the IRAS data alone. The dust 
mass evaluations for the two spirals NGC 4736 and NGC 
4826 and for M31 are in agreement with the average dust 
mass found in the present work, while the gas-to-dust ratio 
is ~ 100, as in the Galaxy. 

The average B-band optical depth along the line-of- 
sight, found by assuming that the total observed dust mass is 
confined within one third of the optical radius, spans a large 
range of values (tb =0.1 to 50, see Fig. 5). The largest val- 
ues of dust mass, optical depth (tb > 1) and of the dust/gas 
mass ratio are shared by objects characterized by a star- 
bursting activity. The inactive "cirrus" -dominated galaxies, 
which are the typical population that bright optical samples 
select, seem to be not much obscured by dust. 
These results are confirmed by an inspection of the actual 
Ab extinction values inferred from the observed ratio of op- 
tical to far-IR bolometric emissions. Relatively high values 
of Ab (> 1) are displayed only by a minority of the star- 
bursting galaxies (6 out of 28 objects). Note, however, that 
for the latter class of sources the quantity Ab is probably 
quite a poor measure of the real extinction affecting the op- 
tical galaxy, since most of the far-IR emission there comes 
from starburst regions occupying very small volumes in the 
galactic body. In these cases too the average extinction is 
probably not larger than 1 mag. 

There are some intriguing effects to notice in Fig. 7. 
The starbursting activity seems to imply a heavier dust en- 
richment of the ISM with respect to inactive spirals. The 
M g /Md median values are ~ 1000 for the latter and 400 for 
the starbursts. The ongoing starburst may well have prepro- 
cessed the ISM in these objects significantly longer or more 
efficiently than inactive galaxies did. 

The somewhat low M g /Md value found for the Galaxy in 
Figure 7 with respect to other cirrus dominated objects in 
our sample may be explained by the fact that the value of 
the ratio for the Galaxy mostly refers to the inner region, 
where most of the dust is concentrated, because such re- 
gion is more strongly weighted in the observations of the 
Galactic emission. On the contrary, observations of exter- 
nal galaxies include the contribution of the outer dust-poor 
environments. So, we probably have to take into account a 
systematic bias when comparing gas-to-dust ratio measure- 
ments for the Milky Way and for distant galaxies. 
Finally, there seems to be a trend favouring higher values 
of the mass in dust Ma relative to the total gas content 
M g in higher luminosity normal inactive spirals. Assuming 
that this trend just reflects an increased metallicity, this ef- 
fect parallels a similar trend of higher metallicities at higher 
masses inferred for galaxies of any morphological types (see 
a recent review in Edmunds and Terlevich 1992). 
A number of abundance indicators, from the specific Mg2 to 
the more general and widespread color-magnitude relation- 
ship, show that larger fractions of supernovae debris are re- 
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tained at greater galaxy masses (probably due to the deeper 
gravitational wall). An effect of this kind may be seen in 
our observations of the dust fraction in IR galaxies. The 
fact that starbursting galaxies in Fig. 7 do not follow the 
same relationship may be due to a much different history of 
metal production. 
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Table 3. Physical parameters of the two-component dust model for the average galaxy spectrum (°) 



cirrus component 
(Xc=Hx©) 



Grain type 



T ^peak 

(K) (H 



mass fraction 
(%) 



starburst component 
(X C =12O X0 ) 



*peak 



(K) (H 



mass fraction 
(%) 
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Notes to Table 3 



(°) The average galaxy has a spectral shape as in Fig. 43, a distance of 70 Mpc, and a flux 5ioo = 5 Jy. The total dust masses in 
the cirrus and starburst components are 2 10 7 M e and 3.4 10 5 M e , respectively. 
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IRAS position 




TV* J. 

Distance 




optical photometry* 


name 




R.A. 




6 




D 


A 

A e 


D 25 


m B 


morp. type 






(h) 









/TV/I"— ~\ 

(Mpcj 


( ) 


(") 


mag 




0028-22 


00 


28 


38.0 


-22 


53 


38 


106.8 


21.1 


67 





14.6 


SB(s)b?pec 


0030-25 


00 


30 


10.8 


-25 


53 


02 


37.6 


21.9 


80 





14.3 


SAO 


0030-22 


00 


30 


48.9 


-22 


38 


08 


245.0 


13.5 


45 





16.3 


SBbc 


0040-23 


00 


40 


17.2 


-23 


49 


57 


83.2 


15.7 


52 





14.4 


SB(r)a 


0045-25 


00 


45 


5.2 


-25 


33 


40 


3.6 


473.0 


1799 





8.2 


SAB(s)c 


0100-22 


01 


00 


22.0 


-22 


38 


09 


471.2 










Sy2 


0116-24 


01 


16 


22.0 


-24 


12 


21 


94.0 


22.9 


61 





14.2 


SB(r)ab 


0128-22 


01 


28 


03.7 


-22 


55 


40 


21.2 


95.5 


234 





11.7 


SAB(rs)c 


0201-23 


02 


01 


35.2 


-23 


33 


14 


65.6 


21.6 


75 





14.1 


(R)SB(r)bc 


0206-23 


02 


06 


59.9 


-23 


39 


04 


70.8 


32.4 


110 





13.2 


Sbc 


0222-25 


02 


22 


49.5 


-25 


00 


53 


45.2 


47.3 


132 





12.5 


SB(s)cd 


0223-24 


02 


23 


42.4 


-24 


56 


03 


138.8 


20.4 


65 





15.2 


(R)SB(r)a 


0300-23 


03 


00 


23.7 


-23 


03 


43 


18.5 


115.0 


313 





11.3 


SB(r)c 


0311-25 


03 


11 


22.0 


-25 


54 


29 


22.4 


94.4 


237 





11.8 


S... 


0318-25 


03 


18 


53.5 


-25 


41 


29 


19.2 


18.2 


65 





13.6 


Sb 


r\ o o r* n a 

0335-24 


03 


35 


19.7 


-24 


39 


48 


19.9 


66.1 


219 





11.6 


SB(s)cd 


0338-22 


03 


38 


45.6 


-22 


43 


29 


20.8 


55.0 


188 





12.5 


(R)SAB(s) 


0346-22 


03 


4b 


47.5 


-22 


16 


59 


56.0 


15.5 


72 





13.8 


(R)SBa 


0408-23 


04 


08 


14.7 


-23 


44 


46 


56.4 


35.9 


129 





13.8 


SO/a 


0427-26 


04 


27 


28.4 


-26 


49 


12 


55.2 


26.9 


77 





13.5 


SB(r)ab 


0433-25 


04 


33 


35.0 


-25 


14 


05 


61.6 


18.6 


81 





14.9 


Sb 


f\ A O l"7 n A 

0437-24 


04 


37 


00.9 


-24 


16 


52 


56.0 


21.4 


67 





14.0 


Sb 


2101-23 


21 


01 


16.5 


-23 


38 


03 


160.0 


10.0 


16 


8 


16.8 


merger? 






40 


44.2 


-25 


35 


UD 


a 7 a 
4 1 .0 


on n 

ZU.U 


100 





14.0 


Sb 


2155-25 


21 


55 


26.5 


-25 


07 


32 


121.1 


20.0 


51 





15.1 


(R)SBa 


2210-26 


22 


10 


18.8 


-26 


23 


49 


65.2 


42.0 


123 





13.2 


SA(s)0/a 


2217-24 


22 


17 


57.9 


-24 


55 


50 


63.2 


40.0 


117 





12.7 


(R)SA(r)0 


2219-23 


22 


19 


35.1 


-23 


20 


13 


107.2 










S r 


2233-26 


22 


33 


00.8 


-26 


18 


31 


18.9 


82.0 


274 





11.5 


SAB(rs)bc 


2235-26 


22 


35 


56.3 


-26 


06 


38 


48.4 


46.0 


173 





12.8 


SA(s)ab 



* optical photometry taken from ESO-LV data base 
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col. (4) and (5) no beam aperture correction 
col. (6) and (7) beam aperture correction with a mm = a a /3 
col. (8) and (9) beam aperture correction with a mm = a a 
•k mm photometry taken from Chini et al. 1984 



